Correlated electrons in optically-tunable quantum dots: 
Building an electron dimer molecule 
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We observe the low-lying excitations of a molecular dimer formed by two electrons in a GaAs 
semiconductor quantum dot in which the number of confined electrons is tuned by optical illumina- 
tion. By employing inelastic light scattering we identify the inter-shell excitations in the one-electron 
regime and the distinct spin and charge modes in the interacting few-body configuration. In the 
case of two electrons a comparison with configuration-interaction calculations allows us to link the 
observed excitations with the breathing mode of the molecular dimer and to determine the singlet- 
triplet energy splitting. 



PACS numbers: 73.21.La, 73.20.Qt, 73.43.Lp, 31.15.vn 

The creation and control of electrons confined in semi- 
conductor quantum dots (QDs) enables fundamental 
studies of Coulomb interaction at the nanoscale and al- 
lows to define optimal architectures for solid-state quan- 
tum information processing pQ. In the case of two elec- 
trons, singlet and triplet spin states with a well-defined 
energy splitting are under active investigation for the im- 
plementation of logic gates for quantum computing [THS] . 

In the dilute limit of few electrons, a paradigm shift 
occurs as the liquid-like properties of the electrons - 
captured by a single-particle picture- cease to dominate 
and the impact of Coulomb correlation favors the for- 
mation of a molecular electron state with strong radial 
and angular order [3H9]. The molecular dimer formed 
by two electrons is especially relevant being the bench- 
mark for the exploration of Coulomb correlation at the 
nanoscale. This state has been the subject of several 
theoretical investigations (e.g. [T0HT3] ) albeit difficult to 
realize experimentally since it requires soft confinement 
potentials that are comparable to Coulomb interaction 
strengths. Two-electron correlated QDs were studied in 
magneto-transport experiments at finite fields [T4] . 

Like in atomic nuclei [15] , the emergence of strong cor- 
relation among the QD electrons manifests in the ener- 
getic structure of the low-lying few-body modes that are 
accessible by inelastic light scattering [9] [16] . Optical in- 
vestigations of the characteristic low-lying excitations of 
two electrons in the QD remain challenging. Such ex- 
periments are enabled by combining spectroscopy tools 
with charge tunability with single electron accuracy, as 
obtained in devices containing self-assembled InAs QDs 
[T7( [T8] where two-electron singlet and triplet states were 
observed by inelastic light scattering [19]. In these self- 
assembled QDs, however, the large confinement energy 
prevents the formation of a molecular dimer state. 

In this Letter we report the identification of a dimer 



molecular state of two electrons at zero magnetic field 
in nanofabricated GaAs/AlGaAs QDs. The evidence is 
found in the low-lying spin and charge excitations ob- 
served by resonant inelastic light scattering. Crucial to 
these experiments is the continuous charge tunability by 
photo-depletion that allows the creation of QDs with a 
selected number of electrons. We first identify the one- 
electron regime characterized by a single monopole exci- 
tation between Fock-Darwin (FD) levels — the states of 
the two-dimensional harmonic oscillator- that are simply 
linked to the in-plane confinement energy. 

The observation of additional distinct spin and charge 
modes as electron numbers in the QDs are increased as 
we sweep the illumination intensity signals the emergence 
of a ground state composed by two and more interact- 
ing electrons. We show that the measured low- lying 
excitation mode energies are incompatible with a the- 
oretical framework that neglects correlation effects. We 
propose an interpretation based on a configuration in- 
teraction (CI) approach (aka exact diagonalization) that 
precisely reproduces the inelastic light scattering spec- 
tra. In the case of two electrons the CI analysis reveals 
links between observed excitations and the breathing and 
center-of-mass oscillations and rigid-rotor rotations of the 
molecular electron dimer. 

In the CI interpretation the rotation energies are 
ft 2 M r 2 el /2/, with / being the momentum of inertia and 
M re \ = 0, ±1, ±2, . . . the angular momentum in the rela- 
tive frame. To comply with the antisymmetry of the total 
two-electron wave function, M re \ must be even (odd) for 
a singlet (triplet) state. Therefore, the energy separation 
between the triplet and singlet ground states is the ro- 
tation quantum ft 2 /2J, that as shown below is accessible 
by the inelastic light scattering experiments. 

The sample is constructed as an array (100 fim x 
100 jam) of identical QDs (10 4 replica to improve the 
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signal-to-noise ratio displaying the same optical prop- 
erties [20] and inter-dot spacing of 1 fim). It is fab- 
ricated by electron beam lithography and inductively 
coupled-plasma reactive ion etching on a 25 nm wide, 
one-side modulation-doped Alo.1Gao.9As/GaAs quantum 
well. The measured low-temperature two-dimensional 
electron density and mobility are 1.1 • 10 11 cm -2 and 
2.7 • 10 6 cm 2 /Vs, respectively. SEM images of such QDs 
are shown in Figs. [TJa) and (b). 

The inelastic light scattering experiments were per- 
formed in a backscattering configuration (q < 2-10 4 cm -1 
where q is the wave- vector transferred into the lateral di- 
mension) at T = 2 K. A tunable ring-etalon Ti: Sapphire 
laser with a wavelength of 789 nm in resonance with the 
QD absorption (whose onset is at 813 nm) was focused 
on the QD array with a 100 jam diameter area, and the 
scattered light was collected into a triple grating spec- 
trometer with CCD detection. 

Photo-depletion of the QD electrons was achieved 
by continuous-wave HeNe laser excitation at 633 nm 
that creates electron-hole pairs in the AlGaAs barrier. 
Whereas the photoexcited electrons in the AlGaAs bar- 
rier contribute to charge compensation of the ionized 
donors, the photoexcited holes are swept in the QDs flat- 
tening the parabolic in-plane potential (i.e. decreasing 
the confinement energy Tiuoq and contributing to the re- 
moval of the confined electrons by electron-hole recombi- 
nation [2lJ[22]) I c f- mse t to Fig. [TJd)] . As shown below 
the method allows to fully photo-deplete the QDs and 
tune their population with single-electron accuracy by 
changing the HeNe laser intensity /HeNe- 

Distinct spin and charge intershell monopole modes 
(with AM = 0, M is the total angular momentum) were 
detected as a function of /HeNe by varying the linear po- 
larization of the scattered photons [parallel (perpendic- 
ular) to that of the incoming laser for charge (spin)]. 
The condition of photo-depleted QDs (number of con- 
fined electrons N = 0) is defined by the absence of any 
inter-shell excitations. This is achieved at /HeNe = 1-82 
W/cm 2 as displayed in Fig.JjJd) [23]. The N = 1 regime 
is characterized by a single excitation at 2huo (degener- 
ate in spin and charge) between FD levels with the same 
angular momentum. Therefore its energy position mea- 
sures directly Tiujq = 1.1 meV, as shown in Fig.JTJc). This 
condition is reached at ^HeNe = 1-46 W/cm 2 . Peculiar to 
this regime is the absence of multiplets at higher energy, 
as confirmed by calculation of the light scattering cross 
section within a perturbative approach under the reso- 
nant conditions employed here [24] . 

Additional spin and charge modes appear at higher 
energy [see Figs.^a) and [4] as we further reduce /HeNe- 
The emergence of multiplets in the spectra experimen- 
tally defines the occurrence of two electrons in the QDs 
at /HeNe = 0.68 W/cm 2 . In this regime of shallow con- 
finement, huo ~ 1-2 meV, the electrons are highly corre- 
lated [6] [3 [9] . To identify the degree of correlation of the 
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FIG. 1: (Color) (a) SEM images of a QD pillar and (b) of 
the QD array, (c) and (d) Resonant inelastic light scattering 
spectra of spin (fitted with red line) and charge (blue line) 
excitations for N = and N = 1. The N = 1 intershell tran- 
sition is shown in the inset to (c). The horizontal segments 
are the FD levels. The sketch of the QD in-plane parabolic 
potential is shown in the inset to (d). The photo-generated 
holes (filled black circles) accumulate in the depletion region 
close to the physical borders of the etched pillar provoking a 
flattening of the confining potential (dotted curve). Removal 
of the electrons (filled light-blue circles) occurs via photo- 
recombination with the photo-excited holes (red arrows). 



two electrons we theoretically evaluate the spectra em- 
ploying the full CI method [25]. This provides accurate 
energies and wave functions of strongly interacting sys- 
tems, by bulding TV-body states as linear combinations 
of many Slater determinants. The latter are obtained by 
filling with N electrons in all possible ways a truncated 
FD level basis set (here truncation includes the lowest 
10 shells, with a relative error on excitation energies of 
« 10 -4 ). From CI wave functions of ground and excited 
TV-body states we evaluate inelastic light scattering ma- 
trix elements and spectra [9] [16] . 

The spectra predicted by CI calculations for N = 2 
[Fig. [2^b)] are remarkably similar to the experimental 
ones. In the evaluations the only free parameter, Tiujq, 
is adjusted to fit the energy of the four peaks shown in 
Fig. [2|a). The spectra appear as two doublets — each 
formed by two peaks in the charge (blue lines) and spin 
(red lines) channels, labeled respectively as Ci, Si and 
C2, S2 in Fig. [5|b). Both spin excitations of Fig. [2|a) 
are slightly red shifted (« 0.2 meV) from their charge 
counterparts. These splittings are also seen in the theo- 
retical calculations in Fig.^b), which also reproduce the 
relative intensity of the peaks in each channel. The value 
of hujo = 1.6 meV used in Fig. [2^b) is larger than that 
measured at TV = 1, as expected for a lower /HeNe [22 . 

The down triangles in Fig. [2^b) summarize instead the 
results of a self-consistent Hartree-Fock (HF) calculation 
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FIG. 2: (Color) (a) Experimental spectra (black lines) of spin 
and charge excitations at IneNe 0.68 W/cm" 2 . Red and 
blue lines are gaussian fits, (b) CI spectra for N = 2 (the 
peaks are artificially broadened by gaussians with FWHM = 
1.5 meV). The down triangles point to the HF modes for 
huoo = 2.154 meV. (c) Comparison between CI (dashed lines) 
and WM modes (solid lines) with houo = 1.6 meV. 



[T3] of spin and charge transition energies that neglects 
correlation effects. This approach was extensively used to 
interpret previous light scattering spectra in the many- 
electron regime [26-28]- By fitting the position of Si 
[lowest down triangle in Fig. [2^b)] to its measured value 
[Fig. [5|a)] by properly tuning hujo, then the remaining 
three HF excitations are blue shifted more than 1 meV 
from their measured values. On the other hand, by fit- 
ting the position of the higher HF doublet (not shown) 
then Si is red shifted by ~ 1 meV. This analysis there- 
fore rules out an interpretation in terms of single-particle 
transitions between FD shells and establishes the onset 
of electron correlation as we increase the electron popu- 
lation from N = 1 to N = 2. 

The impact of Coulomb correlations is represented by 
the dimensionless radius r s of the circle whose area is 
equal to the area per electron (in units of the Bohr ra- 
dius [16 ). For the CI calculation of Fig. |2fb) r s = 3.39. 
In this dilute regime we expect the two electrons to form 
a molecule. This is confirmed in Fig. [2|c) by comparing 
the CI excitations of Fig. [2^b) [dashed lines in Fig. [2|c)] 
with the analytical predictions for a perfectly formed 
"Wigner" molecule [11] (WM) [solid lines in Fig. [2jc)]. 
The WM excitations have a composite nature: Ci (C2) is 
a single (double) excitation of the breathing mode, with 
energy Tiujb (2Tiuub), Si is a dipolar center-of-mass ex- 
citation hujQ plus a rotational quantum ft 2 /2I (the sys- 
tem as a whole keeps M = 0), and S2 adds a breathing 
mode quantum. Remarkably, spin and charge excitations 
scale differently with huo, since ujb = a/3^o whereas 
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FIG. 3: (Color) Pair correlation functions for the CI states as- 
signed in Fig. [2] The unit length is -£qd = (h/m^ujo) 1 ^ 2 . (a) 
Conditional probability P(x, y; xo, yo) of measuring an elec- 
tron in the xy plane provided another one is fixed at position 
(xo,yo — 0), labelled by a black dot, where xq is located at 
the average value of r. From left to right: ground (GS) and 
excited states Ci and C2, respectively. The squares' size is 
7x7, and the 15 equally spaced contour levels go from blue 
(minimum) to red (maximum), (b-c) Pair correlation func- 
tions g(r) vs r for the ground state, charge (b), and spin (c) 
excitations, respectively. 
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The WM estimate for h 2 /2J, and hence the 



singlet-triplet splitting, is 0.517 meV. 

The agreement between CI and WM predictions in 
Fig. |5|c) is very good — with small deviations only for 
the higher excitations — implying that electrons in the 
ground state freeze their relative distance (cf. Fig. [3] 
and Supplementary Discussion). The formation of this 
Wigner molecule, which develops after the emergence of 
rotational bands [9] as r s increases, occurs relatively soon 
for the dimer since it involves only nearest-neighbor par- 
ticles (cf. Supplementary Discussion). 

In Fig. [3]^a) we plot the conditional probability 
P(x, y\ xo, yo) of finding an electron at position (x,y) 
provided the other one is fixed at (xo,yo) (black dot). 
In the ground state (GS, left panel) the well developed 
correlation hole around the fixed electron suggests the 
freezing of the relative distance r. This is best seen from 
the pair correlation function g(r) [black line in Fig.^b)], 
which for N = 2 is the square modulus of the relative- 
motion wave function. Since g(r) is the probability that 
two electrons are at distance r, its well developed max- 
imum around r ~ 1.5 [in units of the characteristic QD 
lenght ^qd = (h/m*^) 1 / 2 ] identifies the inter-electron 
equilibrium distance of the molecule. Besides, the plots 
of Ci and C2 in Fig. [3^ a) [Fig.[3|b)] show respectively one 



4 



charge T=2K 




2 4 6 8 10 12 14 
Energy shift (meV) 

FIG. 4: (Color) Evolution of spin and charge excitations as a 
function of JneNe- Curves are best-fits to the data with zero, 
one or two gaussians. Up triangles point to the CI energies of 
spin (red) and charge (blue) modes. 



and two additional nodal surfaces in the xy plane (nodes 
along r), pointing to the excitation of respectively one 
and two quanta of the breathing mode. Note that the 
molecule is not perfectly formed due to residual electron 
derealization [g(r = 0) =^ in Fig. |3jb)]. Spin exci- 
tations [red lines in Fig. [3^c)] have a different nature: 
g{r — 0) = for symmetry [g{r) ex r Mrel as r —> and 
M re i = 1] and Si (S2) has no node (one node) at finite r, 
consistently with the the excitation of no (one) breathing 
mode quantum. A similar analysis in the center-of-mass 
coordinates confirms the assignments of Fig. [2|c). 

Figure [4] summarizes the evolution of the measured ex- 
citations as a function of IiieNe- The data show that: (i) 
hujo increases at lower IueNe since the modes are blue 
shifted [22]; (ii) N increases for I^eNe < 0.68 W/cm 2 
since the two low-lying spin and charge modes do not 
split significantly, contrary to the expected behavior for 
N = 2 [29]. Indeed the spectra at J H eNe = 0.16 W/cm 2 
agree well with the CI excitation energies for N = 3, 
shown as up triangles in Fig. [4] The best agreement 
with the experimental data is obtained with r s = 2.18 
(houo = 2.8 meV), a value significantly smaller than the 
one for N = 2. This decrease of r s weakens electron cor- 
relations and it stops at iiieNe = 0. Then the spectra re- 
semble those obtained in similar QDs studied previously 
by us [9j [16] and interpreted as due to the low- lying spin 



and charge excitations of N > 3 electrons with r s ~ 2. 

In conclusion, we have reported the characteristic exci- 
tations of two correlated electrons (a molecular dimer) in 
nanofabricated GaAs quantum dots. Charge tunability 
with single electron accuracy is obtained by photodeple- 
tion. These studies open new opportunities for under- 
standing the ground states of few electrons in semicon- 
ductor quantum dots in the regime of strong Coulomb 
interaction and weak confinement potential. 
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Supplementary Discussion: Theoretical analysis of 
Wigner localisation in a broad range of r s 

Since the electron dimer is a finite system, there is 
no clear-cut transition to the WM state as r s increases. 
A good indicator of the crossover to the WM exploits 
the irrelevance of the spin degree of freedom in this clas- 
sical limit [4] [6] . Specifically, the pair correlation func- 
tions g(r) for the ground and Si states shown in Fig. 3(c) 
must overlap in the WM limit. In fact, these represent 
the square moduli of the relative-motion wave functions. 
Since in the WM limit Si and the ground state differ only 
for the excitation of one rotational rigid-body quantum 
plus a center-of-mass oscillation (cf. main text), the rela- 
tive motion — and hence g(r) — is left unchanged. How- 
ever, the two states have different spin symmetries: Si 
and the ground state are a spin triplet and singlet with 
one node and no node for g(r) at r = 0, respectively. 
Therefore the g(r)'s are generically different and overlap 
only in the WM limit, when all the wave function weight 
is removed from the origin. Consistently, electrons fully 
localise in space. 

This crossover to the WM is analyzed quantitatively by 
evaluating the functional distance d(r s ) between the pair 
correlation functions g(r) of the ground state [#gs(0] 
and Si [gs^r)}: 

poo 

d(r s ) = 2ir dr r \g GS (r) - g Sl (r)\ . (1) 
Jo 

d(r 3 ) is plotted as a function of r s in Fig. |5j Note that 
the evolution is smooth [HJ [9], with d(r s ) —> only for 
r s — >• 00. Nevertheless, one may locate the crossover 
to the WM state at r s « 4 where the slope of d(r s ) 
significantly changes. This estimate is consistent with 
the theoretical analysis of Ref. 6, Note the proximity of 
the experimental value (filled black circle in Fig. [5| to 
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of g(r). A feature peculiar of the WM is the parabolic 
dispersion of the rotational band, ~ ft 2 M r 2 el /27, with 
the momentum of inertia / univocally fixed by r s . This 
may be seen in Fig. |6j where we compare the CI lowest- 
energies of the dimer for each M re \ sector (black circles) 
with the parameter-free WM rotational band (continuous 
red line). The agreement is satisfactory up to M re \ = 2 
and then deteriorates as M re \ increases. The rationale is 
that the WM melts as the energy (M re \) increases, the 
larger r s the higher the energy. According to Ref. [lOj the 
results of Fig. [6] show that the WM is "floppy" for this 
value of r s . 



FIG. 5: Functional distance d(r s ) vs r s as obtained by CI 
calculations. The unit length is £qd = (h/m* uq) 1 ^ 2 . The 
continuous line is a guide to the eye. 




FIG. 6: CI lowest-state energy of the electron dimer vs rela- 
tive angular momentum M re i for r s = 3.39. The continuous 
red line is the lowest rotational band of a perfectly formed 
WM (with an arbitrary offset). 

the crossover region, which confirms our interpretation 
in terms of the WM dimer. 

In a previous work [9] we have demonstrated the emer- 
gence of rotational bands for N = 4 at relatively high 
density (r s = 1.71). This transition to a correlated 
molecular state is distinct from the formation of the WM, 
which implies electron localisation and occurs at larger 
values of r s [50] . This may be seen in Fig. [5j where the 
value of d(r s ) for N = 2 corresponding to the same den- 
sity as in Ref. 9 (red filled circle in Fig. |5| is significantly 
distant from the crossover region. 

The assignment of rotational bands in Ref. 9 identi- 
fies states with different values of M re \ and same profiles 
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